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Abstract: We study the quantization of the curved spacetime created by ultrarela- 
tivistic particles at Planckian energies. We consider a minisuperspace model based on the 
classical shock wave metric generated by these particles, and for which the Wheeler - De 
Witt equation is solved exactly. The wave function of the geometry is a Bessel function 
whose argument is the classical action. This allows us to describe not only the semiclassical 
regime (S — ► oo), but also the strong quantum regime (S — > 0). We analyze the interac- 
tion with a scalar field <fi and apply the third quantization formalism to it. The quantum 
gravity effects make the system to evolve from a highly curved semiclassical geometry (a 
gravitational wave metric) into a strongly quantum state represented by a weakly curved 
geometry (essentially flat spacetime). 
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1. Introduction 

The spacetime geometry created by ultrarelativistic sources, that is, the gravitational 
shock wave spacetimes, have arisen great interest recently^ 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ] . These 
backgrounds are relevant to describe the particle scattering at the Planck energy scale. 
Quantum scattering of particle fields and strings by this class of metrics have been studied 
in the approximation in which the geometry is treated classically, i.e., as a background 
field. 

A further step in this direction is to quantize the geometry itself. As it is known, so 
far there is not a quantum theory of gravitation to fully carry out this program. Although 
a conventional quantum field theory of gravitation lacks to be renormalizable, information 
about the quantization of the spacetime geometry can be obtained by solving the Wheeler 
- De Witt equation! 13 ] : 



tf( (3) £)=0. (1) 



In the canonical description of General Relativity, the space-time metric has the 3 + 1 
decomposition 

d 2 s = (N 2 - N.N^dt 2 - 2N l dx l dt + gij dx l dx J , (2) 



that is 



-N~ 2 N~ 2 
( .V -'.V //'•' .V •-'.Y,.Y, ) • (3) 



In the region between two space like hypersurfaces t = ti and t = tf, the Einstein equations 
determine the sequence of three-geometries on the space-like surfaces of constant t, the 
dynamical important object being the 3-geometry &>Q. The dynamics of the gravitational 
field is entirely described by the so called Hamiltonian constraint H. = 0. In the quantum 
theory, this becomes an equation for the state vector = 0, which takes the form 
of the functional differential equation (1). There are also the others constraints of the 

classical theory, but at the quantum level, they express merely the gauge invariance on 

(3) 

Classically, one can know both ^G(to) and Q(to) at some time parameter to and 
determine the 4-geometry ^Q, but quantum mechanically, one can only know ( 3 )(?(t ) 



2 



or G(to) and therefore, one has a certain probability for ^Q(t). The manifold of all 
possible - the so called superspace - in which each point is a metric gij(x), has the 
metric 

Gijki = Vg(9ik9ji + QiWjk - 9ij9ki) , (4) 

with signature (— , +, +, +, +, +). In order to solve Eq. (1), one considers spacetime 
symmetries and restricts the degrees of freedom to be quantized (minisuperspace models). 

2. Pure Gravity Model 

Let us quantize now the gravitational shock wave geometries. The minisuperspace 
metric to quantize is 

ds 2 = —(N 2 - N u N u )dv 2 + 2N u dudv + F(p, u)du 2 + dx 2 + dy 2 (5) 

where u = z — t, v = z + t are null variables and p = \/ x 1 + y 2 . The classical shock wave 
metrics have the form of Eq. (5) wittJ 10 ] 

F(u,p) = f{p)5{u) , (6) 

and the Lagrange multipliers taken N u = 1 , iV = 0. This expression of the metric is 
generic for any ultrarelativistic source; its particularities entering only in the form of the 
function f(p). This expression can be extended to D - dimensions and to shock waves 
superimposed to curved backgrounds I 4 ' 7 !. Eq. (5) can also represent a sourceless plane 
gravitational wave. In such case 

F(x,y,u) = (x 2 -y 2 )D(u) , 

D(u) being an arbitrary function of u. 

The metric on superspace has the components: 



Cam 


- I p3/2 




2 ' 






= G1133 = 


2 ' 


G2222 


= ^3333 = 


— G2233 = 
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the others components vanish. 

The three - curvature is given by 

< 3 >i* = -g^ign + ^(^ 11 )- 2 (V^ 11 ) 2 = ^[-V 2 ± F + ^(V ± F) 2 ] 

where the subscript _L refers to the (2, 3) spatial transverse coordinates. 
Thus, the Eq. (1) in our case reads, 



= . (7) 



Here, for the sake of convenience, we have taken a 3+1 decomposition with respect to the 
hypersurface v = const, that will play the role of time. Notice also that the supermomen- 
tum constraint is already satisfied classically, since for the metric (5), one has R vu = 0. 
Thus, we do not need to consider the quantization of this equation! 14 ]. 

Let us consider now a minisuperspace model by freezing out all the transversal degrees 
of freedom and quantizing the longitudinal one, i.e. we write, 

F(u,p) = f(p)D(u) , (8) 

where f(p) is the profile function characterizing the classical shock wave metrical . The 
function D(u) will represent the degree of freedom to quantize (classically, D(u) = S(u)). 

By discretizing the variable u — > u n and making the change of variable s 2 (u n ) = 
D(u n ) we can transform the functional differential equation (7) into the following ordinary 
differential equation, 



£ + fi p-l).-.iU«7(,) 



*(s„) = , (9) 
where p accounts for the arbitrariness in the operator ordering, and 

C{p) = -2 (vlf - ^ (V±/) 2 ) = 2f WR(p) . (10) 
The whole wave function will read 

9(D) = Y[9(D(un)). 



For p = 1 (which corresponds to the Laplacian ordering prescription), Eq. (9) can be 
brought into a Bessel equation of index v = 0. We choose the solution that remains finite 
when D — > 0. Thus, our solution reads, 

*(F n ) = J (2 v /C( / 9) J D(w n )) = Jo ^2^2 (3)i2F(p,u n )) . (11) 

We can obtain the semiclassical limit by considering large arguments of the Bessel 
function in equation (11). Thus, 

-1/2 



~ ^2 WRF(p,u n )^ cos 2^2 WRF{p,u n ) - tt/4 



(12) 



On the other hand, the semiclassical regime can also be directly studied from the 
Hamilton - Jacobi equation for our system, 

-\ fv2 {§y + Fi/2 (3) ^)=°' 

where S is the Hamilton - Jacobi principal function that by direct integration reads, 

S = ±2\/2 WRF . (13) 

We see that the argument of the wave function (11) is the action of the classical 
Hamilton - Jacobi equation. Thus, for large S, the Bessel function has an oscillatory 
behavior and the wave function Eq. (11), is given by 

* ~ 4 exp{iS/h} + § exp{-iS/fi} . (14) 

that is, Eq. (12) gives the right semiclassical limit. 

Semiclassically, the wave function \1/ is picked around the classical metric geometry. 
It appears natural to interpret the oscillatory behavior Eq. (14) as describing ^ in the 
classical allowed regime. The classically forbidden regime, instead, appears with a real 
exponential behavior (Euclidean signature region). 

For small S, (S h), the behavior is non - oscillatory. In particular, the Bessel 
function reaches a finite value in the limit S — > 0, i.e. Jo(0) = 1. This can be interpreted 



as a genuine quantum behavior, as opposite to the large S sector which describes the 
semiclassical regime. Notice that 5^0, i.e. D(u) = 0, is flat spacetime and that this 
appears in the strong quantum regime. For F > the action is real and Eq. (14) describes 
\& in the classically allowed regime. For each S, i.e. each F, ^ describes a classical 
configuration. All configurations are classically allowed. 

There, | W | 2 (which to some extent can be seen as proportional to a probability), is 
a maximum. We will come back to this interpretation by the end of the paper. 

For > the action is real and Eq. (12) describes the classically allowed regime. If 
( 3 )R < then ty(S — > oo) ~ ex.p{S} / S . This is an exponential growing of the semiclassical 
wave function. It is the opposite of what happens in the case of tunnel effect, and like 
the situation of the falling of a particle into a well potential, indicating the presence of an 
instability. That is, the classically forbidden regime does not corresponds to tunnel effect, 
but to an unstable (exponentially growing) behavior. 

3. Gravity plus matter 

We will include a scalar field in the analysis of the wave function. The Wheeler - De 
Witt equation in this case reads 



here it is understood that the discretization of the variable u is already made in an analo- 
gous way to that of the previous section. 

To make an analysis of this equation we use the following decomposition in modes 
of the scalar field motivated by the solution of the Klein - Gordon equation, i.e. the 
semiclassical regime 



S 2 
5F 2 




5(f) 2 



+ (<P, U f + - [(c/>, p f + m 2 <p 2 } U(F, 0) = . (15) 




where —k + k- + k 2 = — 



o 

m . 
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Thus, replacing jU and 4> }P into Eq. (15) we can write approximately (up to a neg- 
ligible logarithmic term proportional to In | (f>k/4>k |)> 



+ F 2 



5F 2 5(j) 2 
where the effective potential £4(F, (p) is given by 

2 ( 3 )i2(p) 



-t/ fc (F,0)|*(F,0) = O 



(16) 



+ [k 2 _ - 2k+k-F + k\F 2 + (k 2 - m 2 )^- 1 ] (j) 2 , 



(17) 



If we use the following transformation of variables 

X = Fsinh(0) , T = Fcosh(0) , 

this equation can be interpreted as a Klein - Gordon equation for our wave function 
with a time - dependent potential Uk{X,T). When Uk < the wave function \1/ will 
have a real exponential behavior, this corresponds to a classically forbidden regime, in 
which geometries have Euclidean signature. For Uk > 0, instead, the wave function has 
an exponentially oscillatory behavior corresponding to a Lorentzian signature (classically 
allowed) regime. 

The analysis of the potential can be done directly in terms of the variables F and (j). 
To this end we define: 



(j) 2 k 3 



2 (3) R+ ( k 2 -m 2 ) 



k+F 
~k~7 



Thus, 



The zeros of the potential can thus be obtained from the condition 



(18) 



"o = A)(i - A)) 2 , 



(19) 



and the extrema of the potential from 
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(20) 



Let us consider now the potential as a function of F for slices of <fi = constant. There 
are two main cases: 

a) a < 

In this case the potential Uk will always be bigger than zero having a minimum at 
a value (3 m > 1 (solving the third order equation (20)). Thus, we have here the normal 
oscillatory behavior of the wave function. 

b) a > 

Here we have three subcases depending on the value of a: 

i) a > 8/27: The potential grows monotony from — oo at f3 = to +oo at (3 — > oo 
passing through zero at (3q > 1. 

ii) 4/27 < a < 8/27: In this case the potential shows a local maximum and minimum, 
both in the region Uk < 0. For (3 > (3q > 1, Uk > 0, and we have a steady growing. The 
maximum and minimum are located at < (3 max < 2/3 and 2/3 < (3 m i n < 1 respectively, 
and they can be obtained by from Eq. (20). 

iii) < a < 4/27: This case presents two regions where Uk > 0. These two regions 
are Pi < (3 < /3$, where < $ < 1/3 and 1/3 < /3$ < 1; and (3 > (3% > 1. The two local 
maxima and minima are of course located at (3q < (3 max < (3q and (3q < (3 m i n < (3q. 

The dependence of Uk on <p is clearly a growing parabola (see Eq. (17)). The minimum 
value is reached at the origin </> = and it will eventually cross a zero of the potential iff 

( 3 )# < 0. 

Thus, we have seen that we can classify the two regimes: Euclidean (exponential 
behavior of the wave function) and Lorentzian (oscillatory behavior) according to the sign 
of the effective potential Uk- A weak potential Uk corresponds to the semiclassical regime. 
This is precisely what happens when F is large, i.e. S — > oo, namely a strong curvature 
regime, for instance a state picked around a singular gravitational wave. A strong potential 
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Uf~ corresponds to a truly quantum regime in which 5 ^ 0, i.e. F — > 0, that is, a smooth 
and weak curvature state. 

If we prepare the system in a classical configuration, for instance \<f>\ = 1 and F — > oo 
(that would correspond to the semiclassical problem studied in Ref. 6,10 of a Klein - 
Gordon field in a classical shock wave geometry), we observe that we fall into the case 
a), i.e. Uk > 0, since ^'R > 0. Due to the particle creation (that we will consider in 
the next section) of the field 0, the system evolves towards the minimum of the potential 
Uk making F small and |0| big, that is, the system evolves towards a quantum regime of 
smooth and weak curvature. 

4. Third Quantization 

In this section we study the possibility of third quantizing the wave function of 
the ultrarelativistic particles. A motivation for the third quantization is to overcome the 
problem of negative probabilities just as was the case for the Klein - Gordon equation (see 
for example Ref. 15). 

Let us then consider the Eq. (15) for an ultrarelativistic particle and an homogeneous 
massless scalar field 0, that up to the operator ordering ambiguity reads, 

dt -dl + 72 WRexp{-6t} \ 0) = , (21) 

where we have made the change of variables (at each point u n ) 

t n = -iln(F n ), (22) 

which allows the variable t n to run from — oo to +oo and to give to it a time - like 
interpretation. 

The general solution to Eq. (21) can be written as 

®j(t, 0) = II (V8 Wife 3 *") exp{zj</>} , (23) 

n 

where Z v is a Bessel function of first or second kind. 



It is interesting to note here that the case of pure gravity is recovered for j = 0. As 
j is related to the energy of the matter field 0, it seems that the presence of matter fields 
excite the gravitational modes and allow for particle production as we shall immediately 
see. 

Following Ref. 16 we can define in states proportional to the Bessel function 
J2ij (^/S ( 3 )i?e 3 *™^ for t — > — oo (which are natural positive frequency modes in the in re- 
gion), and out states proportional to the Hankel function H^j (^/S ( 3 )i?e 3t ™j for t — > +oo. 
Then, one can compute the Bogoliubov transformation coefficients between these two ba- 
sis, and since the Hankel function is a linear combination of positive and negative subindex 
(frequency) Bessel functions one obtains particle production of the outgoing modes with re- 
spect to the in vacuum. It appears that the spectrum of produced particles has a Planckian 
distribution at temperature, T = l8^8nG/3, of the order of the Planck temperature. 

Thus, the interpretation of our results can be made in terms of the creation of ultra- 
relativistic particles carrying with them its own geometry. The interaction with the matter 
fields deplete the gravitational energy bringing the ingoing semiclassical state, for instance 
picked around a shock wave metric, (that is, a strongly curved geometry) into a state in 
the quantum regime represented by a weakly curved geometry (essentially flat spacetime) . 
The evolution of the system is from the classical into the quantum regime. The initial 
semiclassical configuration is a highly curved geometry, the final configuration is a weakly 
curved one. 
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